Localness of energy cascade in hydrodynamic 
turbulence, II. Sharp spectral filter 



We investigate the scale- locality of subgrid-scale (SGS) energy flux and inter- 
band energy transfers defined by the sharp spectral filter. We show by rigorous 
bounds, physical arguments and numerical simulations that the spectral SGS flux 
is dominated by local triadic interactions in an extended turbulent inertial-range. 
Inter-band energy transfers are also shown to be dominated by local triads if the 
spectral bands have constant width on a logarithmic scale. We disprove in particu- 
lar an alternative picture of "local transfer by nonlocal triads," with the advecting 
wavenumber mode at the energy peak. Although such triads have the largest trans- 
fer rates of all individual wavenumber triads, we show rigorously that, due to their 
restricted number, they make an asymptotically negligible contribution to energy 
flux and log-banded energy transfers at high wavenumbers in the inertial-range. We 
show that it is only the aggregate effect of a geometrically increasing number of 
local wavenumber triads which can sustain an energy cascade to small scales. Fur- 
thermore, non-local triads are argued to contribute even less to the space-average 
energy flux than is implied by our rigorous bounds, because of additional cancel- 
lations from scale-decorrelation effects. We can thus recover the -4/3 scaling of 
nonlocal contributions to spectral energy flux predicted by Kraichnan's ALHDIA 
and TFM closures. We support our results with numerical data from a 512'^ pseu- 
dospectral simulation of isotropic turbulence with phase-shift dealiasing. We also 
discuss a rigorous counterexample of Eyink (1994), which showed that non-local 
wavenumber triads may dominate in the sharp spectral flux (but not in the SGS 
energy flux for graded filters). We show that this mathematical counterexample 
fails to satisfy reasonable physical requirements for a turbulent velocity field, which 
are employed in our proof of scale-locality. We conclude that the sharp spectral 
filter has a firm theoretical basis for use in large-eddy simulation (LES) modeling 
of turbulent flows. 
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Abstract 



I Introduction 



This paper is the second part of a study of the scale-locahty properties of the energy 
cascade of three-dimensional turbulence. In the first part'^ (hereafter referred to as I) 
we investigated energy transfer defined by means of low-pass and band-pass filters for 
smooth, dilated kernels, which permit simultaneous resolution of the physical processes 
both in space and in scale. We demonstrated there the scale-locality of interactions 
involved in turbulent energy transfer, by rigorous analysis and by numerical simulation. 
See alsol^. The discussion in these papers did not apply to sharp-spectral filters, which 
do not satisfy the modest decay conditions in physical space that were employed there. 
The origin of scale-locality with the sharp spectral filter is, in fact, a rather more subtle 
issue. The problem is of importance, however, since most numerical studies of turbulent 
scale-locality employ sharp spectral filters. We believe that misinterpretation of the 
numerical results has led to a number of misunderstandings in the literature. 

Early papers on turbulent energy cascade which employed spectral analysis — such as 
those of Obukhovl^, Onsagerl^, and Heisenber^^' — argued for scale-locality. Those works 
proposed a "cascade process" in which two modes of similar wavenumber transfer energy 
to a mode with twice that wavenumber, implying steps in spectral space of increasing size 
with increasing wavenumber. This picture was later confirmed by detailed closure calcu- 
lations of Kraichnan, both for his (Eulerian) direct-interaction approximation (DIA)'^, his 
abridged Lagrangian- history direct-interaction approximation (ALHDIA) and the test- 
field model (TFM)'^. In the late 80's and early 90's, however, this traditional view was 
contested in numerical simulation studies by a number of authors: Brasseur & Corrsin'^, 
Domaradzki & Rogallo'^, Yeung & Brasseurf^l^ and Ohkitani & Kida'^. Those works 
presented results suggesting that energy cascade is dominated by highly nonlocal triadic 
interactions, with one very energetic mode at the lowest wavenumber ko catalyzing trans- 
fer between two modes at high wavenumbers k < k', for k' — k = 0{ko). This transfer 
process was described as "local transfer through nonlocal interactions" . In contrast to the 
traditional picture, the energy flow through spectral space was suggested to proceed via 
steps of a fixed, small size /cq, driven by low- wavenumber advection. Chasnov^S! argued 
that such nonlocal, advective sweeping interactions should be incorporated into spectral 
large-eddy simulation (LES) via a stochastic force that produces random backscatter of 
energy. If such views were correct, they would present a fundamental challenge to the 
Kolmogorov picture of local energy cascade and the universality of small-scale turbulence. 
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Several theoretical papers published shortly thereafter defended the traditional view 
of turbulent scale-locality™SIM2l2l] particular, WaleS(^, Zhou^Sl and EyinkEH 
all argued that scale-locality is recovered when Fourier modes are suitably combined or 
averaged. Zhou'^^'^ showed numerically that locality does hold, when some necessary 
summations are made, and he verified the quantitative prediction of Kraichnan'^^'^S ^j^at 
non-local contributions to energy flux decay as a power oc s^^'^ of the "scale-disparity 
parameter" s = A^min/^mod (where A;inax, ^med, ^min are the maximum, median, and mini- 
mum wavenumbers of a triad, respectively). It seems, however, that this conclusion was 
subjected to doubt in a subsequent investigation by Zhou et al. (1996)'^ which indicated 
that the effects of the largest scales are significant. The paper of Eyink ( 1994)^1 proved 
rigorously that spectral energy flux, if averaged in wavenumber over an octave band, is 
dominated by local triads. His argument yielded a rigorous upper bound 0(s^/^) on the 
nonlocal triadic contributions, which is larger than Kraichnan's scaling prediction but 
still vanishing as s ^ 0. However, the paper of Eyink also proved that the spectral en- 
ergy flux, without additional averaging, may be non-local. He constructed a velocity field 
with Holder continuity properties analogous to those observed in turbulent flow for which 
the instantaneous spectral energy flux is dominated by nonlocal advective interactions, 
somewhat similar to the effects seen in the numerical simulation^ ^^ * ^^ * ^ ^^^. 

The debate about the locality of interactions contributing to energy flux with the 
sharp spectral filter has continued unabated. Recently, the issue was addressed by nu- 
merical simulations at much higher Reynolds numbers by Alexakis et al. (2005)^^ and 
Mininni et al. (2006,2008)'25M1^ ^ho concluded, in agreement with the earlier simu- 
lations ^ * ^^ * ^^ * ^^ ', that the important interactions contributing to interband transfers are 
those among highly non-local wavenumber triads with one leg at the energy peak. On the 
other hand, they concluded that the spectral energy flux, involving a suitable summation 
of Fourier modes, is dominated by local triads. Domaradzki & Carati (2007,2009)'232SI 
confirmed these results and furthermore numerically calculated Kraichnan'^^^'221 locality 
function W{s), a quantity which measures the fractional contribution of nonlocal triads. 
They verified Kraichnan's prediction that W{s) oc s^/^ for s <^ 1, in agreement with ear- 
lier studies of Zhou'^^'^. On the other hand, they found no difference between the sharp 
spectral filter and various graded filters, in apparent contradiction with the analytical 
work of Eyink'21'2^ who drew a clear distinction between the sharp spectral filter and 
smooth, dilated filters. 

It may be an understatement to say that no coherent picture of the scale-locality of 
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turbulent energy cascade is presented by the current literature on the subject. Open 
issues that continue to be discussed include: the need for additional averaging (over 
space,time, ensembles, etc.), qualitative distinctions between sharp spectral versus smoothly 
graded filters, the quantitative fraction of energy flux contributed by non-local triads, 
and possible differences between spectral binning with logarithmic versus linear bands. 
Even the proper definition of "energy flux" is unclear, since the definition used in the 
rigorous proofs with graded filters'^^^'^ is the "SGS flux" familiar from large-eddy simu- 
lation modeling, which involves additional advective subtractions to define the subscale 
stress which are not present in the conventional definition of the spectral energy flux. 

The present paper has two main purposes. Our first goal is to present new results 
on scale-locality of spectral transfer quantities, including rigorous estimates, physical 
arguments and numerical simulation data. In particular, we shall show rigorously that 
both energy flux and inter-band transfers defined by the sharp spectral filter — without 
any additional averaging — are dominated by local triadic interactions. Our proof of 
locality rests on four ingredients: (1) The SGS flux (defined in I and in the next section) 
as the unique Galilean invariant measure of the inter-scale energy transfer at each point 
in the flow, (2) scaling properties that are observed empirically to hold for the turbulent 
velocity field, (3) wavenumber conservation, which constrains the number of Fourier 
modes contributing to the energy flux, and (4) the essential use of "logarithmic" spectral 
bands whose width increases proportional to wavenumber. The second aim of this paper 
is to attempt to collect the disparate results in the literature into a clear and consistent 
picture of the local cascade process in three-dimensional turbulence. We shall try to 
provide answers to many of the open issues mentioned above which continue to be debated 
in the literature. 

To aid this latter goal, it may be useful to summarize here the essential reason for 
the dominance of local triadic interactions in energy cascade. It is true, as indicated 
by numerical studies^ ^ * ^^ * ^^ * ^^ * ^"^ * ^^*^^ ', that individual non-local triads make a much bigger 
contribution to energy flux than individual local ones, due simply to their having one of 
their modes at the energy containing scales. However, such comparisons of single triads 
vastly underestimate the aggregate contribution from local triads. Whereas the number 
of highly elongated, nonlocal triads grows moderately with increasing wavenumber, the 
number of local triads increases much more rapidly. The net contribution to energy flux 
from the exploding population of local triads dwarfs the contribution from the smaller 
number of nonlocal triads. While this paper was being written, we discovered a work 
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of Verma et al)^ which arrives at essentially the same conclusion, by means of a non- 
rigorous perturbative closure analysis of energy transfer. They stated there that "...the 
shell-to-shell energy transfer rate is found to be local and forward. This result is due to 
the fact that the nonlocal triads occupy much less Fourier space volume than the local 
ones." A related observation was also made by Alexakis et al.^^, who observed that the 
fractional contribution of non-local triads to energy flux is reduced "since many more 
local triads contribute in the global summation." In fact, these ideas, as we shall discuss, 
are already implicit in the arguments advanced by Kraichnan for scale-locality, using 
his ALHDIA and TFM closures'^^'^. Our work supports these conclusions and extends 
them, by rigorously bounding the nonlocal contributions to spectral flux, without any 
heuristic approximations. Our analysis implies as well the dominance of local triads to 
inter-band energy transfers, for wavenumber bands of constant width on a logarithmic 
scale. We also explain by physical arguments how decorrelation effects (as discussed in I 
for the graded filter) further diminish the nonlocal contributions, giving agreement with 
Kraichnan's asymptotic scaling predictions. 

We verify our theoretical analysis by analyzing the velocity field u generated from a 
direct numerical simulation of the incompressible Navier Stokes equation 

+ (u ■ V)u = - Vp + z/V^u + f (1) 

V-u = 

which is solved pseudo-spectrally in a periodic box of 512^ grid-points. Here, p is the 
pressure, u is the viscosity, and f is an external stirring force. We advance in time using 
the second-order Adam-Bashforth scheme and employ the phase-shift method to remove 
aliasing effects'^. The fluid is stirred using Taylor-Green forcing: 

f = /o[sin(fcjx) cos{kfy) cos{kfz)-k — cos(fc/x) sm{kfy) cos{kfz)y] (2) 

where /o = 2 is the force amplitude, and kf = 2. The viscosity u is 0.87 x 10^'^ and the 
Reynolds number Re = UrmsL/v ^ 11800, where L = 27r is the size of the simulation 
box and u^ms = y/'^Etot = 1-63 is the rms velocity. The Reynold's number based on the 
Taylor scale A = 2ny/E^/ ( J dkk'^E{k)^ is Re\ = 615. We analyze a time snapshot 
of the flow after it has reached a statistically steady state. The energy spectrum E{k) of 
the flow, shown in Figure [l| has a reasonable k~^^^ scaling until around k = 40, before 
dissipation effects start to dominate. The mean energy flux = {^k) is shown in 
Figure |2| 
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Figure 1: Energy spectrum for the simulation being analyzed. The forcing is at kf — 2. 
A straight hne with slope —5/3 is added for reference. 




Figure 2: Average energy flux (n^) from the simulation. 
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The outline of our paper is as follows. Our main new results are presented in the 
following Section II. Subsection II. A contains some preliminary definitions and discussion. 
In subsection II. B we establish ultraviolet locality of spectral SGS flux, and in subsection 
II. C we establish infrared locality. In the next section III, we examine how our work 
relates to previous research on scale-locality, in particular the analysis of Kraichnan'^ 
using spectral closure and recent numerical studies^ ^ * ^^ * ^^ * ^^ * ^^ *^. Section IV formulates 
our final conclusions. Two appendices contain some technical material. Appendix A on 
our rigorous proofs and Appendix B on Kraichnan's ALHDIA closure. 

II Main Analytical and Numerical Results 
A Preliminaries 

In a periodic domain = [0, 2tt)^, one may represent a velocity field u by its Fourier 
series expansion: 

u(x) = J2 u(k)e^''-. 

Restricting this sum to wavenumbers satisfying |k| < K yields the low-pass filtered field 
u*~^, while restricting the sum to |k| > K gives the complementary high-pass filtered 
field u^^ = u — u<^. We shall also employ in our discussion band-pass filtered fields 

for K < Q. We employ a special notation [K] for dyadic (octave) bands [K/2,K] and 
denote the corresponding band-pass velocity field by u^^l The reader may assume here 
that |k| denotes the Euclidean norm |k| = \/k1 + k'^ + k'^. In fact, due to very sub- 
tle aspects of harmonic analysis, our proofs do not apply with full rigor for the Eu- 
clidean norm, although they are rigorous for the norms |k|i = + -|- j/c^l or 
|k|oo = maxjlfc^:!, \ ky\, \ kz\}. See Appendix A for a careful discussion of these dehcate 
mathematical issues. In fact, we believe that our arguments do apply with the Euclidean 
norm on wavenumbers and that the abstract mathematical counterexamples that com- 
promise their rigor have different properties from those enjoyed by the turbulent velocity 
fields. In support of this belief, all of the numerical results presented later shall employ 
the Euclidean norm and, as may be seen, they agree with our theoretical analysis. 

Just as with any other filtering scheme, one may employ the sharp spectral filter to 
decompose the Navier-Stokes equation ([T]) into large-scale and small-scale components 
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and then consider the corresponding energy balances. Following a standard analysis, 
presented in I, one arrives at the sub-grid scale (SGS) energy flux 

n^(u, V, w) = -d.uf' (^{v,w,)<'' - vf^'wf^^ (3) 

as the representation of the transfer rate of energy from modes with wavenumbers < K 
to modes with wavenumbers > i^' at each space point. We have distinguished here 
a straining mode u, an advecting mode v, and an advected mode w, although, in the 
physical balance equation, u = v = w. This permits us to formulate the property of 
scale-locality in a precise manner. As in I, we say that energy flux is ultraviolet (UV) 
scale-local if replacing any single u or any combination of u's in 11/^ (u, u, u) by u^^l gives 
an asymptotically negligible contribution for P ^ K. Likewise, we say that energy flux is 
infrared (IR) scale-local if replacing any single u or any combination of u's in ni^(u, u, u) 
by u[*^l gives an asymptotically negligible contribution for Q <^ K. 

Most of the previous studies of scale-locality with the sharp spectral fllter have em- 
ployed what we call an "unsubtracted flux" 

n^- = -(9,M,)<^M,) or - (9,M,)<^M,)<^. 

These two quantities have the same average over space as the SGS energy flux, because 
the sweeping terms which distinguish them integrate to zero: f dx diuj^ uf^ u"^^ = 0. 
However, the pointwise values are different and the locality properties are different as 
well. Note that the quantities nx(u, u, ul'^]) and n^"^(u, u, ut^l) = -diuf^{uiuf^)^^ 
are not the same even after space-averaging, since the sweeping effects no longer integrate 
to zero: J dx. diU"^^ uf^ {v^^^)^^ ^ 0, as was recognized by Domaradzki and Caratil^. 
Indeed, it is easy to see that the quantity 11^** cannot be scale-local instantaneously, at 
each point in the space domain. If we boost the flow with a uniform velocity uq (a A; = 
mode), then 11^* will develop a contribution directly proportional to this velocity. Since 
Uq can be made arbitrarily large, triads involving the A; = mode can give the dominant 
contribution at any flxed K. This argument does not apply to the space (or ensemble) 
average of 11^'^, which is Galilei invariant. As we shall discuss later, the mean value 
is in fact scale-local, as concluded by early workers ^ * ^ * ''' * ^ * ^ * '^" *, although its locality 
properties are less robust than those of the SGS flux Ili^. We examine only the latter for 



all of our detailed estimates in this section, and we return to discuss 11^'^ in section III 
Scale-locality of interactions is not a general property of solutions of the Navier- 
Stokes equation, e.g. the nonlinear interactions in laminar flow are dominated by the 
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most energetic modes. Rather, the locahty of the energy cascade in a turbulent flow 
depends crucially on some scaling properties of the solution. The operative scaling laws 
for sharp-spectral filtered quantities are closely related to those for velocity increments 
6u{£; x) = u(x + £) — u(x). Heuristically, 

^>K ^ ^[K] ^ Su{iK), VU<^ ~ 6u{iK)/iK 

with ix = 27t/K the length-scale corresponding to wavenumber K. More precisely, the 
scaling laws that are used in our arguments below are for volume-averages: 



{\u>''\y/P, (|u[^]f)i/P ~ {const.)urUKL)-''', {\Vu<''\y/P ~ (const.) ^(irL)!--^- 

(4) 

with < (Xp < 1 and p > 1, for inertial-range wavenumbers K satisfying KL ^ 1. 
Here L is the integral length-scale, Urms is the root-mean-square velocity, and (.) is a 
volume average. Note that ap = (p/p, where (p is the scaling exponent of the pth-order 
absolute velocity structure function: {[6u{i)]^) ~ i'^p. The ultimate source of these scaling 
properties is empirical evidence from experiments and numerical simulation j^^ * ^"^ * ^ ^. The 
precise mathematical forms of the scaling laws ^ that are used in our proofs are discussed 
in Appendix A. Here we emphasize just one salient point: the scaling law stated for 
(|u['^l 1^)^/^ in Q is valid only for logarithmic bands [K] and does not hold for a band 
[K, K + A] with fixed spectral width A. This point is essential for the proof of infrared 
locality, as we shall see. 

B Ultraviolet Locality 

We shall begin by treating the case of ultraviolet locality, which is somewhat easier and 
allows us to introduce the basic method of argument. 

1 Theoretical Analysis 

We must consider the quantities Hi^(u['^l, u, u), nA'(u, u^^l, u), and Hx(u, u, u'-'^l) for 
P ^ K to establish UV locality in the straining, advecting, and advected modes, re- 
spectively. The first of these is trivial, since di{uP)<^ = ior P > 2K and thus such 
modes give a vanishing contribution to the strain. The second two UV locality proper- 
ties are less trivial but their physical origin lies in the simple fact that modes at smaller 
scales contain less energy. Thus, the contribution from modes P ^ i^' to the stress at 
wavenumber K 

TKin,u) = (uu)<^ - u<^u<^ 
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is small. We will now show this through precise estimates. 

We first consider nx(u, u, u'^^) = —diuf^TK{ui,u^p). When P > 2K the subtrac- 
tion in the stress vanishes, so that there is no difference between the SGS flux and the 
"unsubtracted" flux. This quantity then reduces to: 

The condition that the advecting mode lie in the band [P/2 — K,P + K] arises from the 
constraint that the wavenumbers of the two velocity modes in the stress must sum to a 
value < K. We thus see that, for P ^ K, the two stress modes must lie in essentially the 
same high-wavenumber band, since [P/2 — K,P + K] is the same as [P] = [P/2, P] with 
a little extra padding of size K. This argument apphes with equal force to the quantity 

Thus, we may consider these two quantities together. 

A rigorous estimate which establishes UV locality of the SGS flux follows from the 
Holder inequality: 

(|9,<^(t.l^-^'^+^^t.fl)<^|> < (const.) (|Vu<^|^>^/='(|u[f-^'^+^l|^>^/='(|u[^]|3>^/l 

See Appendix A for details. Together with the scaling laws ^ this gives 

(|ni,(u,u,u[^l)|) < (const.) ^(Lir)i-'^='M2^,(LP)-2'^^ ~ £(const.)(Li^)^-=^'^« {K/Pf'K 

(5) 

where e = u^^^/L is the energy injection rate. For any scaling exponent a-^ > 0, the 
factor {K/ PY"^ becomes very small when P ^ K, implying that the high P modes make 
asymptotically little contribution to the space-averaged energy flux. In fact, it is known 
empirically that (J3 = 1/3, the K41 value. Thus, we obtain a bound close to 0{{K / P)"^/^). 

One difficulty in the above estimation is the factor [LKY'^"^ . Since 0-3 < 1/3, this 
factor grows slowly with increasing K and has the potential to render our bound useless 
at very high K, deep in the inertial range. This would be the case if the bound were larger 
than the total energy fiux e itself, so that it would cease to provide a tight constraint 
on the large-P contribution. Since the growth in K is so slow, however, it can be 
easily compensated by taking P large enough. More precisely, define the wavenumber 
P^{K) = i^(/s:L)(i-3'^«)/2a3 > Tj^gj^ upper bound (|5|) is < e for P > P^{K), 

implying that such modes make an asymptotically negligible contribution to the mean 
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flux. This proves rigorously the UV-locality of the space-average energy flux. Note, 
however, that the bound ([s]) is probably far from optimal, as discussed more below. 
The same argument as above can be applied to prove that 

for any p > 1 (Appendix A). In the limit p — > oo, the quantity (|ni4'(u, u, ut^])^)^/^ 
increases to sup^ |ni4:(u, u, u^''^])!, the maximum over the domain, and ap decreases to 
/imin, the minimum Holder exponent of the velocity field. Therefore, by taking larger p, 
our rigorous estimates hold more uniformly over space but also imply less rapid decay 
with increasing P. These results are similar to those established by Eyinkl^SEl and in I for 
graded filters, but weaker. In particular, we do not have pointwise estimates of energy 
flux at each space point x in terms of the local Holder exponent, using spectral filters. 

Up until this point, our discussion has been mathematically rigorous. It is often 
the case, however, that rigorous proofs do not yield the optimal results on complex 
physical problems. Our upper bound 0{{K / PY"'"^) on the mean flux contributed by 
high-wavenumber modes P ^ K is larger than the asymptotic scaling result oc (K/Py^^ 
predicted by Kraichnan ^^° l ^^l Just as discussed in I for graded filters, we can argue on 
heuristic, physical grounds that the smaller contribution found by Kraichnan is due to 
cancellations that arise in the average over space, which are neglected in our crude upper 
bound. For details, see I and*^. The result is that 

(H/^(u,u,ul^]) ~£:(const.)( 

If we assume also the K41 value = |, then we recover the 4/3 scaling prediction of 
Kraichnan'^^'^ who obtained both the fractional contribution {^p/^kY^'^ and the K41 
value C2 = I from his ALHDIA and TFM closures. Both of these results also agree with 
the (^p/^k)^"^^ estimate made by L'vov and Falkovichl^, if one assumes C2 = |- In 
principle, however, all of these results are slightly different due to intermittency effects 
and experimentally distinguishable. 

2 Numerical Results 

We shall now present results from our 512^ DNS, both to investigate the sharpness of 
our rigorous bounds and to test the validity of the physical arguments. 
We first present simulation data for the following quantities: 
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All of our rigorous estimates apply to these two objects, but neither can experience the 
additional cancellations invoked in our heuristic argument, because of the added absolute 
values. Notice that we have included absolute values even inside the low-pass filter < K. 
As discussed in paper I, substantial cancellations could otherwise occur in the space- 
integral which defines that low-pass filter. We have calculated space-averages of the 
quantities (|6| for K = A and for logarithmic bands [P] = [P/2,P] varied continuously 
over values from P = 8 to P = 240, where k^ax = 241 is the maximum wavenumber in 
the simulated flow field. We plot the results in Figures [3^, b, respectively, normalized by 
(lllifl). The two quantities are very similar and both show a slightly faster decay than 
our prediction. This is consistent with the findings of Domaradzki & Carati'^^'^, who 
also see faster decay of UV nonlocal contributions than predicted by Kraichnan'^^'^. This 
is plausibly attributed to the extreme shortness of the inertial range and contamination 
by viscous effects. The decay rate beyond P = 40 is definitely increased by viscosity. We 
have fitted power-laws over the range P = 8 to P = 40, which is the putative "inertial- 
range" of our DNS based on the scaling of the energy spectrum (Figure [T]). Over this 
range we obtain a p~o-83 scaling for the first quantity in ^ and P~^-^^ for the second, 
both distinctly faster than P~^/^. 

We next present results for the mean quantities nif(u, u, u^^') and nx(u, u^^l, u), 
again for K = A and for P = 8 — 240, normalized by the mean flux Hk = These 
data are plotted in Figure |4] and show good agreement with our predictions. When 
fitted over the "inertial range" from P = 8 to P = 40, we obtain a P~^-^'^ scaling for 
ni<'(u, u, ut^]) and p^^-^^ ni^(u, u'"^', u). This is a slightly faster decay than P~^^^, 
again in accord with the recent results by Domaradzki & Carati'^^SHI. Note that the plots 
for nj^(u, u, u'-^l) and n/^(u, u^^l, u) are identical for large P. This makes good sense, 
because, as mentioned above, the maximum and middle wavenumbers can differ at most 
by K and thus those two modes arise from nearly the same wavenumber range. Another 
interesting feature is that for P > 110, the partial flux becomes negative. This shows up 
as a kink on the log- log plot, a feature that was also observed by Domaradzki & Carati'^. 

C Infrared Locality 

We turn now to a discussion of the infrared locality of the spectral energy flux. We shall 
prove, by similar arguments as in the previous section, that the large-scale components 
of all three velocity modes — the straining, advecting and advected modes — contribute 
negligibly to spectral energy flux. The importance of the large-scale advecting mode con- 
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Figure 3: For K = A, we plot (a) {\diuf^ ' ^uf^\) \)/{\nK\) and (b) 

{\diuf^{\uf^uf~^'^^^\y^\)/{\UK\) using logarithmic bands [P/2,P]. The straight 
lines with —2/3 slope are for reference and extend over the fitting region which gives 
slopes of -0.83 and -0.98 for (a) and (b), respectively. The two plots are essentially 
identical for large P, as expected. 
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Figure 4: For K = A, we plot the global quantities (a) nii'(u, u, ut^^)/!!^ and (b) 
n^(u, ut^], u)/!!^^ using logarithmic bands [f,-P]. The straight lines with —4/3 slope 
are for reference and extend over the fitting region, which yields slopes of -1.54 and -1.42 
for (a) and (b) respectively. The two plots are essentially identical for large P. The kink 
at P ~ 110 in the log-log plot is due to a change in sign at that wavenumber. 
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tribution has, in particular, been a major source of contention in the hterature. Several 
studies^ ^ * ^^ * ^"^ * ^^ * ^^ * ^^ * ^^ * ^^ ' employing numerical simulations have concluded that nonlocal 
advective interactions are primarily responsible for the energy transfer. In this picture, 
energy cascades to smaller scales by very many small steps in wavenumber, of fixed size, 
mediated by the largest, energy containing modes. If true, such a picture would have 
dramatic implications for Kolmogorov's concept of small-scale universality and on the 



theory and practice of LES modeling. However, we argue below (and in section III ) that 
this picture is false, and is an artefact of the use of constant (unit) width spectral bands 
to define energy flux and energy transfer. 

1 Theoretical Analysis 

We derive exact bounds on the quantities Uxi^^'^K u, u), nx(u, u''^], u), and Hxi^, u, ut'^l) 
for Q K, thus establishing IR locality in each velocity mode. 
We first consider the contribution of the straining mode: 

n^(u[«l, U, U) = -d,uf^TK{u,, Uj). (7) 

The origin of IR-locality here is fairly obvious, because the strain from low-wavenumbers 
is weak. This fact may be expressed as a rigorous upper bound by using the Holder 
inequality: 

\d,uf^TK{u,,u,)n < ( IVuIQir ) ( |r;,(u,u) 



with ^ = ^ + ^- Because TK{ui,Uj) ~ 6ui{iK)Suj{£K), it is expected that 

l/s 



TK{u,u)\'j ~ {const.y,^,{KL)-P^ (8) 

with ps = 2a 2s- See Meneveau & O'Neil'^ and our own figure [6| where we verify this 
scaling for a sharp spectral filter with s = 2. If this result is combined with the scaling 
relation (|4]) for the velocity-gradient, we obtain 

IP \ ^ ( i- \^""'^ ( T r>\^-'^r 2 



\d,urW{u,,u,)Y J < (const.)^(Lg)^-- ui^.iKL)-''^ 

= (const. )e(Q//^)l-"'■(f^L)l-"'■~''^ (9) 

As long as a,. < 1, the factor (Q/i^)^"'^' is decaying for Q K, implying IR locality 
in the straining mode. Just as in the discussion of UV locality, larger p corresponds to 
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estimates more uniform in space. In the IR case, however, the decay estimates improve 
for increasing p because ap is non-increasing in p (Frischl^, section 8.4). For example, 
with r = s = 2p, we see that the exponent 1 — increases with p and the bound on the 
nonlocal contribution becomes tighter. 

Also as in the UV case, there is an additional overall factor (i^'L)^^'^''"''^ which causes 
our bound to deteriorate for KL ^ 1. For example, for p = 1, r = s = 2, (T2+P2 = 0"2+2cr4 
and the latter is slightly less than 1 by the concavity of (p as a function of p. (See, for 
example, Frischl^S, section 8.4). However, our bound ^ is still useful (less than e) if we 
take Q small enough to offset the slight growth in K. Defining the wavenumber Q^, (K) = 
(l/L)(KL)2'^4/(i-a2) g^pj^ ^i^g^^ ^ Q^{K) < K, one need only take Q < Q^{K) to 

ensure that the contribution of straining modes at wavenumber Q make a negligible 
contribution to energy flux. 

We now discuss the large-scale contributions of the other two velocity modes which 
appear through the turbulent stress, i.e. the advecting and advected modes. These two 
cases are far more delicate. We shall consider in detail the case of the advecting mode, 
which has been the source of most of the controversy in the literature. Thus, we shall 
examine 

(10) 



ni^(u,u[^i,u 



for Q < K/2. The case of the advected mode is quite similar and will be discussed in 
parallel. The key to the whole analysis is the following fundamental identity: 



nK(u,u[«U) = -5^«; 



<K 



Ua Ua 



'IV. 



To derive this relation, note that any advected modes Uj with wavenumbers < K — Q m 



(llOl) when multiplied by ■uf^'^ will contain only wavenumbers < K. Hence, the low-pass 



filter < K on the first term in the square bracket [■] can be dropped and the contribution 
of all such modes cancels exactly with the similar contribution from the second term in 
the bracket. The additional restriction to wavenumbers < K + Q for Uj in the first term 
arises from the fact that modes with wavenumbers > K + Q when multiplied by can 
only give wavenumbers > K and these modes do not survive the low-pass filter < K. 



The identity (11) expresses a basic restriction on the range of interactions allowed to 



very nonlocal triads, with advecting modes at low wavenumber Q interacting only with 
advected modes in a band of wavenumbers [K ~ Q,K + Q] at most, with width 2Q. A 
similar identity can be derived as well for Hi^(u, u, u^*^!), by the same argument. This 
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restriction on the range of the nonlocal triadic interactions is the origin of the IR locality 
of the SGS spectral flux in the advecting and advected modes. 

We now wish to demonstrate this by a rigorous estimate. Unfortunately, the scaling 
relations Q that we have employed until now are not applicable to estimate the size 
of the band-pass field u^^~^'^~^''^\ As we mentioned earlier and discuss in more detail 
in Appendix A, the scaling of m^^^ in Q follows from the empirically-known scaling of 
velocity-structure functions, but only for [K] = [K/2,K], a logarithmic band. Here, the 
band [K — Q,K + Q] has a fixed width 2Q, independent of K. However, we may appeal 
to another empirical observation on the scaling of the turbulent energy spectrum. A 
great many experiments and simulations have shown that 

Eik,t)r^Cul^^LikLp'^^^^ (12) 

with (2 = 2/3, for wavenumber binning of unit size. That is, the Kolmogorov power-law 
spectrum is observed not just with octave bands but with bands of essentially infinitesimal 
width. This is, in fact, one of the most robust empirical observations on small-scale 
turbulence. From it we infer that for Q <^ K 

pK+Q 

|^[x-Q,x+Q]|2^ = 2 / dkE{k,t) ~ {const.)ul^,{Q/K){KL)~^''\ 
Jk-q 

The important point here is the factor {Q/K) which becomes smaller with decreasing Q. 
To exploit this estimate, we use the 4-4-2 Holder inequality to derive the bound 

H;,(u,u['31,u)|\ < (const.) / |Vu<^|A7 |u[«l|A7 |u[^-Q'^+«f ^ ' 



~ (cOnst.)^(Lii')l-'^* ■ UrmsiLQ)"'' ■ Urms{Q/Kf\KL)--- 

~ (const.)£(LK)i-'^2-2<x4(g/^)i/2-<x4_ (13) 

See Appendix A for more details. Since (J4 is known empirically to be a bit smaller 
than 1/3, the exponent 1/2 — (14 > 1/6 and thus the factor {Q / Ky/"^'"^ becomes small 
for Q <^ K. This proves the IR locality in the advecting mode of spectral SGS flux. 
The result is due to two competing factors. On the one hand, the low- wavenumber 
velocity modes have greater magnitude, reflected in the factor Urms{LQ)^'^^ which grows 
for decreasing Q. On the other hand, the modes with low wavenumber Q are restricted to 
interact with a smaller set of modes at wavenumber K, reflected in the decreasing factor 
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{Q/Kf^. The latter factor wins the competition, implying IR locality. This is one of 
the most important results established in our paper. 

As with our earlier bounds, there is a troublesome factor (Li^')^~'^2-2o-4 h^qx grows 
with K. However, defining Q^.{K) = K{KLf''^l^^-^''^'^~^ so that l/L < Q^.{K) < K, we 
see that our bound becomes <^ e for Q ^ Q^{K). Thus, IR locality follows, although the 
result is presumably not optimal. Eyinkl^ (and I) used graded filters to derived a sharper 
bound of the form 0{{Q / K)"^^^), whereas the decay proved here for the sharp spectral 
filter is only like {Q/KY^^. As we shall discuss below, there is some evidence from our 
numerics to suggest that the 2/3 decay rate holds also for absolute spectral fiux and even 
faster decay is expected for mean spectral fiux due to decorrelation effects. 

First, however, let us discuss briefiy how the bound (13) is reconciled with the 1994 
counterexample of Eyink'^H which showed that spectral fiux may be dominated by nonlo- 
cal advective sweeping interactions. The counterexample cannot satisfy the bound (13), 
which would imply instead that the local triads dominate. There is no contradiction, 
however. Eyink's counterexample is a Fourier- Weierstrass-type function given by a la- 
cunary Fourier series, with just two wavenumber modes in each octave band [2^, 2^+^] 
for = 1,2,3, .... It has scaling properties similar to those of a turbulent velocity field 



(like (|4])), but it fails to satisfy the strong condition (12) on the energy spectrum be- 
cause its wavenumber modes do not densely populate Fourier space. In this respect, the 
counterexample appears very different from typical turbulent velocity fields. Thus, it is 
true as a general mathematical fact that energy fiux defined by the sharp spectral filter 
may be dominated by nonlocal triads (unlike fiux defined with graded filters). However, 
this is not the case for the turbulent fields satisfying condition (12). In particular, the 
counterexample of Eyink is only slightly related to the "local transfer by nonlocal triads" 
observed in turbulent simulations^ ^^ * ^^ * ^^ * ^^ * ^^ * ^^ * ^ -^*^. As we shall discuss in section HI 
those DNS observations are completely compatible with locality of energy cascade. 

All of our results on IR locality have, to this point, been mathematically rigorous 
upper bounds. However, we expect much faster decay of non-local contributions for 
the space-average fiux than the bounds proved above, due to cancellation of fiuctuating 
positive and negative parts. It is possible to provide more physical arguments to explain 
what we believe is the true scaling of IR non-local contributions to mean spectral fiux. 
For details, see I and'^. Here we note just the final results that 

_ _ /f \i+Ci 

H^(u[«U,u), Hx(u,u[«l,u)~^(^^j • (14) 
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The result is different for locality in the advected mode, so we give a few details. By 
wavenumber conservation with Q < K: 



The limited number of modes contributing to leads to a major cancellation. 

Indeed, expressing the triple correlation heuristically as increments and using fusion 
rules as in I, we obtain the scaling result ioi Q <^ K 



K I <-Q 



K-Q \ I f-K 



~ e — 



HQ 

C- /„ N 2 + Cl 



e 



The cancellation leads to an additional factor of Ik/^q^ which gives a faster decay for 
this non-local contribution. 

It is worth noting that a similar cancellation appears in the partial flux nx(u, u^'^\ u). 



where it is crucial to give the scaling ~ {d-x / (^qY^^^ in (14). The physics behind IR 
locality in the advecting mode is that when Q is very small, then, due to wavenumber 
conservation, only a few wavenumber modes within distance Q oi K contribute to the 
flux, implying insignificant transfer from those triads. In section [TTTl we shall rederive this 
result within the ALHDIA closure of Kraichnan'^, which embodies in a quantitative form 
the same basic ideas employed heuristically above. As for the UV case, our prediction 
of the IR non- locality correction scaling as {Q/KY'^'^^ can in principle be distinguished 
from the scaling prediction {Q/KY~^'^ of Kraichnanl^ and L'vov & FalkovichliZl, due to 
small effects of intermittency. 

2 Numerical Results 

We now present simulation results on IR locality. We first consider space-averages of the 
absolute values of the partial fluxes. These should be devoid of any cancellation effects 
but are constrained by our rigorous upper bounds. In Figjs] we plot (|nif(u['^l, u, u) |) 
(o), (|n;^(u,u,u[«l)l) (x) and (|ni^(u, ut^l, u)|) (+), for K = 100 and for Q ranging 
continuously over all wavenumbers Q < K/2. Fitting with power-laws over the "inertial- 
range" from Q = 8 to Q = 40 gives scalings of Q^'^'', Q^'^^ and respectively. The 

first is in good agreement with our rigorous estimate {Q/ KY~^'^ , for a2 = 1/3. The 
other two quantities have scaling exponents exactly between the value 2/3 and the 1/6 
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Figure 5: For K = 100, we plot (|n^(u[«l, u, u)|)/(|n^H) (°); (|n^(u, u, u[Ql)l)/(|n;,|) 
(x); and (|nx(u, u''^], u)|)/(|nx|) (+) using logarithmic bands [Q/2,Q]. The straight 
lines are for reference and have a 2/3 slope. They extend over the fitting region which 
yields slopes of 0.77, 0.42, and 0.45 for (o), (x), and (+), respectively. Plots (x) and 
(+) are almost identical, as expected. 



exponent of our rigorous upper bound. In Figj6] we also test the scaling prediction (|8]) 
of the spectral SGS stress for s = 2 and obtain good agreement with the K41 value 
exponent p2 = 2(T4 = 2/3. 

In Fig. |7|we plot -(n;^(u[«], u, u)) (o), (n^(u, u, ul^l)) (x), and (H^^ (u, ul^l, u)) 
(+). Fitting with a power-law over the range Q = 8 to Q = 40, we obtain scaling 
laws of Q^ *^^, Q^'^*^ 1 and Q^'^^, respectively, in fairly good agreement with our heuristic 
estimates Q^^^, Q"^^^ and Q'^^^. The IR decay rates are a little slower than those predicted, 
in agreement with the findings of Domaradzki & Carati'^^'^. This can be attributed to 
the extreme shortness of the inertial range and the relatively smooth velocities at large 
scales in the simulation. Note that (n;f(u[*^l, u, u)) < over the range of the graph, 
becoming positive for Q > 60. This can be understood in terms of detailed balance of 
energy in which the band of wavenumber modes [Q] is receiving energy from modes > K, 
but the amount getting progressively smaller for Q ^ K and more negligible compared 
to the positive total mean flux {Hk)- 
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Figure 6: The rms stress (^Tx{u,u)^^^^/Etotai normalized by the energy of the system 
and plotted as a function of K. The straight line with -2/3 slope is for reference and 
extends over the fitting region which gives a slope of -0.71. 

Ill Comparison with Previous Work 

Our entire analysis of the preceding section considered the SGS spectral energy flux g. 
It is more traditional to base discussions of wavenumber locality in turbulence on the 
triplet transfer function 

T(k,q,p) = Re{iki{u*{k)ui{q)uj{p))}5u,p+q- 

This function can be interpreted as giving the mean rate of energy transfer into wavenum- 
ber mode k from mode p induced by mode q and satisfles the "detailed conservation" 
property'^ 

T(k,q,p)+T(k, -q,p) = 0. 

This transfer function is often integrated over spherical shells to define a quantity T{k, q, p) 
which depends only upon wavenumber magnitudes k,p, q, restricted to values which can 
be assumed by the side-lengths of a closed triangle. In numerical simulations'^^'^^'^lESl 
this quantity is represented by 

TiK, Q, P) = 9,^^-^'^]«p-^'«l«f -^'^1, (15) 
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Figure?: For = 100, we plot -(nx(u[«], u, u))/(n^) (o); (nK(u, u, u[«l))/(ni^) (x); 
and (ni^(u, u))/(nii-) (+) using logarithmic bands [Q/2^Q]. The straight lines are 
for reference and have 4/3 and 7/3 slopes. They extend over the fitting region which 
gives slopes of 1.08, 1.90, and 1.27 for (o), (x), and (+) respectively. 
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using band-pass velocity fields with bands of unit width. This function is defined 

pointwise in space, although usually its spatial average is taken over the flow domain. 
From the triplet transfer T{k,q,p) one can obtain the mean spectral flux as 

i*co roo 

= — / dk dq dpT{k,q,p), (16) 
Jo Jo Jo 

or, pointwise in space, 

nr = - Yl T{K\Q',P') = -d,u<''{u,u,). (17) 

K'<K,Q',P' 

As discussed earlier, this quantity has the same space-average value as the SGS energy 
flux (|3|. However, its pointwise properties are quite different, since it is not Galilei- 



invariant unless averaged over space. Thus, the conventional spectral flux (17), and also 
its variant form —diU^^{uiUj)'^^, are pointwise scale-nonlocal objects, dominated by 
large-scale advection effects. All such sweeping interactions, however, can be associated 
to space-transport of energy, with the SGS energy flux remaining as the unique pointwise 
Galilei-invariant measure of energy transfer to small-scales. See Eyinkf^^ and I. It was 
for this reason that we chose to analyze the spectral SGS flux in the previous section, 
since it enjoys much better scale-locality properties than the conventional spectral flux. 

On the other hand, the conventional spectral flux when averaged over space (or, 
assuming ergodicity, over time or over ensembles) is also dominated by local triadic 
interactions. This was realized by the early pioneers of the subject ^^ * ^*^*^ ' and analyzed 
in detail by Kraichnan'^'ISEl using closure approximations. A careful examination of 
Kraichnan's works shows that his arguments for scale-locality are based on the same basic 
ingredients as our rigorous proofs in the previous section. We shall briefly review here 
the asymptotic analysis in section 3 of Kraichnan (1966)1^ using the ALHDIA closure, in 
order to discuss its relation to our arguments and also to facilitate later comparison with 
the numerical studies^ ^ * -^^*-^^ * -^^*^^*^^*^^*^^*^ ^*^. In order to make direct contact between our 
proof and Kraichnan's analysis it is convenient to consider the quantity 

pK pQ poo 

ni^(u,u[«U) = - / dk dq dpT{k,q,p), (18) 
^0 Jq/2 Jk 

which measures the mean energy flux due to advection by modes in the wavenumber 
band [Q]. Kraichnanf^^ considered a slightly different set of quantities — in his notations, 
T^q{K) and II^^q^K) — or the transfer and flux due to advection by all modes with 
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wavenumbers < Q. He obtained asymptotic expressions for these in the hmit Q <^ K. 



However, his analysis and results carry over essentially to the flux quantity (18). 

An important fact used by Kraichnan in his analysis is that, due to wavenumber 
conservation, this partial flux may be written as 

rQ j-K rk+q 

TIk{u,u^^\u) = - dq dk dpT{k,q,p). (19) 

Jq/2 JK-q Jk 

Introducing dimensionless variables u^vhy k = K — qv, p = K+qu, this may be rewritten 



as 

•1 rl-v 



TIk{u,u^'^\u) = - dqq^ dv du T{K - qv, q, K + qu). (20) 
Jq/2 Jo Jo 

The result (20) is exact, involving no approximation. It shows that the restriction on 
allowed wavenumbers interactions provides an additional small factor of for small q, 
which proves essential to obtain IR spectral locality. 

To proceed further requires an explicit expression for the transfer function. Kraich- 
nan's ALHDIA closure yields the formula 

T{k,q,p,t)=BkpqAkpg [ ds[G{k;t\s)U{p;t\s) - G{p;t\s)U{k;t\s)]U{q;t\s); (21) 

J to 

see equation (2.6) of Kraichnan (1966)^°. Here U{k;t\s) and G{k;t\s) are 2-time La- 
grangian velocity-correlation and mean-response functions. The factor Akpq is equal to 1 
if the wavenumbers k,p, q can be the sides of a closed triangle and is equal to otherwise. 
Finally, the factor B^pg = 4:Tr'^k'^p'^q{xy + z^), where x = cosa,y = cos/?, 2; = cos 7 are the 
cosines of the angles a, /5, 7 that are opposite to sides of length k,p, q in the wavenumber 
triangle. Substituting the steady-state, inertial-range scaling forms 

G{k; t\s) = g{e^'^e'\t - s)), U{k; t\s) = ^e^/^k-^^lh{e^l^k^'\t - s)\ (22) 

one obtains the asymptotic expression for k = K — qv, p = K + qu, and q <^ K 

T{k, q, p) ~ -ew{l - w^)C^IK~*/''q~^/\ (23) 

where w = u + v and / is an integral over the scaling functions g{T),r{T). See Kraichnan'^ 
and Appendix B for details. We just note here one signficant feature of the calculation, 
which is the near cancellation between the two input and output terms in (21) for k ^ 
p ^ K. The physics of this was discussed by Kraichnan'^, p. 1733: 
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"These terms separately give contributions proportional to the energy in the 
wavenumbers < g, not to the mean-square vorticity. ...The input and output 
contributions from low wavenumbers are proportional to energy because they 
represent convection as well as straining. Convection of high-wavenumber 
structures by strongly excited low-wavenumber velocity components implies 
a rapid exchange of phase of the high-wavenumber Fourier amplitudes.... This 



exchange is represented in (2.6) [our (21)] by the large, cancelling input and 



output contributions. The net contribution T<g(fc,t) represents the effect of 
straining alone." 

These remarks are consistent with our observation that the spectral energy flux II^'* in 



(17) is pointwise non-local and dominated by convective sweeping, with such effects only 
cancelling in the average over space. The cancellation between input and output terms 
ioi k p K also yields an additional small factor of q crucial to the final result. 



Substituting the closure approximation (23) into the expression (20) yields finally 
that 

n,,(u,u[«],u)^C£f , g<i^, (24) 

for C = Y^C^I (l — ^^) = 0.503. For details, see Appendix B. We thus recover Kraichnan's 
4/3 scaling law for decay of IR- nonlocal contributions. The factor Q^^^ arises from 
integration of the product q^^^ = q^q^^^^ over the octave band [Q/2^Q]. Just as in our 
rigorous proof in the preceding section, there is a competition between the decaying factor 
q^ which arises from the restriction on wavenumber interactions and the growing factor 
which arises from the increase in energy of low-wavenumber modes. For a similar 
calculation and analysis, see Verma et al.'^, section 3. 
We have focused so far on energy flux but many 

^i^^iQ^^^mmmmm instead 

consider energy transfer between wavenumber bands. Our results also imply scale-locality 



of suitable transfers. Note, for example, that the partial flux in (19) for Q < K/2 equals 



n^(u,u[«i,u) = n^(u[^],uf'^',u'^i) 

i-K j-Q j-P _ 

dk I dq I dpT[k,q,p), 

K/2 Jq/2 Jp/2 

^ -T{[K],[Q],[P]), (25) 

with P = 2K. In the last line we have defined the mean triplet transfer function with 
octave bands [K] = [K/2, K]. Because this quantity equals the partial flux, our previous 
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bounds and scaling results all apply, showing that this quantity for Q <^ -ft' is negligible 
compared with the total transfer rate from adjacent band interactions. 

It is very important for the closure calculation presented in this section and also for 
the rigorous proof in the preceding section that such logarithmic bands be employed. 
When demonstrating IR locality in the advecting mode, we have used the fact that 
modes with wavenumbers Q <^ K interact with a very restricted subset of the modes 
contained in u^^\ resulting in an overall weak contribution to the flux. On the other 
hand, if linear bands are used, the modes contributing to u^^~^'^^ are already restricted 
and, thus, taking Q smaller does not impose any further restriction. Thus, the space- 
average of T{K, Q, P) in (15) must increase for smaller Q simply because u^'^^ gets bigger 



at the larger scales, where most of the energy resides. This can be verified directly from 



the asymptotic formula (23) of the ALHDIA closure for T{k,q,p), which is essentially 
the same quantity. However, T{K, Q, P) and other such objects which compare single 
triads vastly underestimate the contribution from the local triads, which, when taken 
into account cumulatively, swamp the non-local interactions through sheer number. 

Such remarks apply to studies which have claimed to verify a cascade mechanism of 
"local transfer by nonlocal triads" in numerical simulations. For example, several such 
^Q].]^g[i3[ i5 | 24 | 25 | j^g^yg calculated the partially-summed transfer 



T{K,P) = J2nK,Q,P) = ^jd'x d. 



\K-1,K] [P-1,P1 



This can be interpreted as the mean rate of energy transfer into wavenumber band 
[K — 1,K] from band [P — 1, P], mediated by all other wavenumbers. The cited nu- 
merical studies have found that the energy transfer to the i^'th band has largest positive 
contribution from P = K — kf and largest negative contribution from P = K + kf, medi- 
ated by modes at the forcing wavenumber kf. This result was suggested to support the 
dominance of nonlocal triads in energy cascade* ^^ * ^^ * ^^ * ^^l In fact, the result should have 
been expected and is completely consistent with scale-local interactions! Notice that 



this behavior is predicted, for example, by the asymptotic formula (23) from Kraich- 
nan's ALHDIA closure, since its integral over q is clearly dominated by the peak of the 
q-5/3 spectrum at q = kj. However, Kraichnan correctly deduced from this expression 
that cascade proceeds by the local triads, not the non-local triads. This can be seen by 
considering instead the similar quantity 



Tm,[p]) = ^ I d'xd.ufwp, 
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employing octave bands. Our previous discussion applies to this object, showing that 
the non-local triads make a negligible contribution. 

To underline this conclusion, we shall exhibit illustrative numerical results from a 
512'^ simulation using hyperviscosit}^ 

dtu + {u-V)u = -Vp + uioiV^fu + i (26) 

We employed the same forcing ^ as for our viscous simulations but with kj = 1, and 
we took hyperviscous coefficient uiq = 2 x 10~^^. These choices were made to achieve 
a constant mean energy flux Uk = £ over as long a range of wavenumbers as possible. 
See Figure [s^. Just below, in Figure [sjo, we plot T{K,P) for linear bins magnified by a 
factor of 10 (dashed line) and T{[K], [P]) for octave bands (solid line). Both quantities 
are plotted as functions of P for K = 50 and have been normalized by the mean fiux e. 

The results for unit-width bands agree with those of the earlier studie^ ^ * ^^ * ^^ * ^^ !, while 
those for octave bands are similar to those presented by Verma et al.'^. Figure 5 for 
logarithmic bands (base-v^ rather than our base-2). However, it is quite instructive to 
plot the two quantities together. The result for T{K, P) illustrates the contribution of 
"local transfer by nonlocal triads", with the familiar feature of peaks separated by 2kf. 
Notice, however, that the magnitude is miniscule compared both with the mean flux e 
and also with the summed effect of all the local triads, which dominates in T([K], [P]). 
This kind of transfer will become an even smaller fraction of the total for larger K, 



decreasing as ~ K~^^^ according to the Kraichnan estimate (24). The plot of T([ir], [P]) 
shows that transfer into [K] has largest positive contribution from the octave band [P] 
with P ^ K/2 and largest negative contribution from the band with P 2K. This 
demonstrates the multiplicative, self-similar nature of the local cascade. 

Our results rigorously disprove a number of suggestions in the recent literature. Alex- 
akis et al. (2005)^^^ and Mininni et al. (2006,2008)'25E6!have claimed on the basis of 1024^ 
and 2048'^ DNS that the ratio nx(u, u''^^, u)/IIk with Q •::=:. kp asymptotes to a value of 
about 0.10-0.20 for large K in the inertial range. A long plateau within the inertial range 



is ruled out by our rigorous estimate (13) which shows that this quantity is bounded at 



least as 0{£{KL) °') for large with a = (T2 + cr4 — 1/2 = 1/6. According to the estimate 



(24) from Kraichnan's ALHDIA closure theory, the decay rate is actually ~ e{KL)^'^^^, 
which is even faster. Whatever is the explanation of the plateau observed in the reported 
simulation^^lEMB] cannot be an inertial range effect. Indeed, the asymptote observed 
in their numerical data occurs at high wavenumbers K beyond the constant flux inertial 
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Figure 8: Using data from the hyperviscous simulation, we plot (a) the energy flux 
(Hk)- In (b), we plot, for K = 50, as a function of P: (sohd) {diu\I^^UiU^p) / (Uk) using 
octave bands [K/2,K], [P/2,P], and (dashed) 10 x {diuf~^'^^Uiuf'^'^^)/{nK) using 
linear bands. We multiply the latter by 10 for comparison. The dashed-line plot shows 
energy transfer in incremental steps with the distance between peaks being 2, which is 
twice the forcing wavenumber, kf = l.This kind of transfer represents a negligible fraction 
of the flux (Hx), as the plot shows, which will decrease even more (~ K~'^^^) for larger 
K. The solid-line plot has peaks at P = 32 and 88, showing transfer in multiplicative 
steps. The solid-line decays like P^-^^ ^^gj, ^^le range [1,30], quite close to the expected 
power P^/^. There is no similar discernible power-law decay in the ultraviolet limit of 
large P due to dissipation. 
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range and within the pronounced "bottleneck" of their energy spectrum. A Reynolds- 
number dependence of the plateau level has furthermore been reported in the most recent 
paper'^(see Fig. 4 there) 

n^,Ku'M,„) ^ ^=0.6-0.7. 

where Kx is the Taylor-scale wavenumber. Since L/A ~ Re^^'^, this is completely consis- 
tent with Kraichnan's estimate (24), which would yield a scaling law Re~'^/^ for the above 
ratio. Inertial-range theory cannot explain, however, why the fraction stays constant for 
K > Kx. We conjecture that this is an effect of the bottleneck. In Figure |9] we plot 
the same ratio using data from our hyperviscous simulation. It shows a decay very close 
to through the inertial range, until it asymptotes at 2-3% over the region of the 

bottleneck, which is known to be more prominent for hyperviscosity'^^HD^ Note that we 
do not observe a similar plateau in our simulation with regular viscosity, which does not 
have a conspicuous bottleneck (see Figjl]). 

Alexakis et al. (2005)'23! and Mininni et al. (2006,2008)'2SE6l have further speculated 
that their numerical results — which they interpreted to show the dominance of non- 
local interactions to mean transfer — may be due to intermittency and long-lived vortical 
structures. However, the bounds on space-average energy flux in the current work are 
little affected by intermittency, since they involve only low-order scaling exponents cr^ 
with j9 = 2 ~ 4. Our estimates for large p and the pointwise estimates for graded 
filters^ ^ will indeed be modified by intermittency, which has the effect of improving the 
IR estimates and degrading the UV estimates. UV locality will nevertheless hold at every 
space point with a positive Holder exponent of the velocity field. 



IV Conclusions 

We have proved by rigorous estimates that spectral energy flux in three-dimensional tur- 
bulence is dominated by local wavenumber triads. No additional averaging over scales, 
as in Eyink (1994)'^, was required. Furthermore, we showed that spectral locality holds 
without additional cancellations resulting from space averaging, which, when included, 
further decrease the net non-local contribution. Our proof exploited four main ingredi- 
ents: (1) The SGS spectral energy flux, which is the unique Galilei-invariant measure of 
the energy flux across scales. It has the sweeping effects subtracted, which is essential to 
prove infrared locality. (2) Scaling properties that are observed empirically to hold for 
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Figure 9: Using data from the hyperviscous simulation, we plot {IIk{u,u^^,u))/{IIk) 
(o) as a function of K. The straight line with slope —4/3 is for reference and extends over 
the fitting range of the graph, which gives a decay rate of K~^'^'^. Beyond this range, the 
plot asymptotes at 2-3%. The plot of the energy spectrum (solid line) suggests that this 
trend is related to the bottleneck in the crossover between the inertial and dissipation 
ranges of the spectrum. 



30 



the turbulent velocity field. In particular, we employ the strong scaling condition eq.(12) 
on the energy spectrum, which rules out the mathematical counterexample of EyinlJ^. 

(3) Wavenumber conservation, which restricts the range of interactions in Fourier space. 

(4) Logarithmic bands of wavenumber modes which, unlike linear bands, represent the 
cumulative effect of the growing number of local triads. 

We were thus able to explain the results of previous research claiming the dominance of 
non-local interactions in the energy cascade* ^^ * ^^ * ^^ * ^^ l Since these studies used linear bands 
in decomposing the velocity field, they were effectively comparing triads of single Fourier 
modes with each other. The most dominant of these triads are the ones with a large-scale 
advecting mode. However, such non-local triads are few, an idea implicit in the work 
of Kraichnan (1966)^*^ and emphasized more recently by Verma et al. (2005)'^. Thus, 
their effect becomes more negligible as the energy cascades to smaller scales, whereas the 
number of local triads grows geometrically, swamping the non-local effects. Therefore, the 
Kolmogorov picture of local cascade is consistent with the presence of non-local triadic 
interactions, since these become vanishingly insignificant in a longer inertial-range. 

We find that the scale- locality properties of the SGS fiux with a sharp spectral filter 
are close to those with graded filters established by Eyinkl^^and in paper I. Our numeri- 
cal results for the space-average of the SGS spectral fiux are essentially the same as those 
for graded filters, in agreement with the earlier findings of Domaradzki & Carati'^^EZl^ xhe 
locality properties of the SGS spectral fiux Uk are possibly not quite as good pointwise 
as those of SGS fiux with graded filters. We could prove some nearly uniform estimates 
on UV and IR locality, but with decay rates that are probably sub-optimal and with ad- 
ditional ii'-dependent factors. More seriously, we could not prove pointwise estimates on 
the SGS spectral fiux involving the local Holder exponent at the space point x. However, 
we would like to moderate the claim of Eyinld^ that spectral energy fiux "is an inappro- 
priate measure of energy transfer." This is true for the unsuhtracted fiux H^'^, which is 
pointwise non-Galilei-invariant and becomes scale-local only when averaged over space. 
On the other hand, the SGS spectral fiux Tlx is scale-local, even in the absolute sense, 
without cancellations due to space-averaging. In particular, the sharp spectral filter has 
a sound theoretical basis for use in large-eddy simulation (LES) modeling of turbulent 
energy cascade. 

The width of wavenumber bands is more important for scale-locality than the grad- 
ing of the filter kernel. As we have seen, transfer functions for individual wavenumber 
triads or transfer functions with wavenumber bands of fixed width on a linear scale may 
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have the greatest magnitude for non-local triads. The situation is completely different 
for transfer functions with bands of fixed width on a logarithmic scale. The transfers for 
such logarithmic bands are dominated by local triad interactions, which, by their vastly 
greater numbers, overwhelm the contribution of the non-local triads. Only the collective 
contribution of the local triads is large enough to explain the observed flux of energy to 
small scales. In addition, it is only quantities summed over logarithmic bands which can 
be simultaneously localized both in scale and in space. This is a basic fact of harmonic 
analysis, equivalent to the Heisenberg uncertainty principle^ZI. To constitute an "eddy" of 
half the size requires twice as many wavenumber modes to achieve the necessary localiza- 
tion in space. If one believes, as we do, that the essential dynamics of turbulent cascade 
involves modes resolved both in scale and in space, then summation over logarithmic 
wavenumber bands is necessary to its physical description. 

Fourier series, without a doubt, are often a useful tool to analyze turbulent cascades. 
However, individual Fourier modes have no physical significance and it is dangerous to 
think of the elementary interactions as those involving individual wavenumber triads. A 
single Fourier mode is a wave which extends periodically throughout all space, whereas 
incompressible fluids (when constant-density, non-rotating, etc.) exhibit no linear wave 
phenomena. Consider the case of turbulent flow in a wall-bounded domain, like air in a 
room or water in a cup. The velocity field for such a flow may be decomposed into Fourier 
modes by imbedding the flow domain in a periodic box. A single triad then represents 
interacting periodic waves, exchanging energy beyond the solid walls of the domain! This 
clearly has no physical reality. Fourier modes are only one convenient basis, among many, 
that may be used in the scale decomposition of the velocity field, and the physics cannot 
depend on the expansion coefficients of any particular basis. Only the sum over the basis 
elements can have intrinsic meaning. 

We have proved in this work that, asymptotically for long inertial ranges, the energy 
cascade proceeds neither in large strides between widely separated Fourier bands nor in 
small incremental wavenumber steps. Rather, it proceeds in multiplicative steps between 
Fourier bands which are adjacent to each other on a logarithmic scale. However, the 
bounds on the non-local triad contributions to energy flux decay only as power-laws, like 
the 4/3-power law of Kraichnan. The cascade process can thus be accurately described 
as "diffuse"'^ or "leaky"H2l. Very large Reynold's numbers are needed, larger than what 
today's numerical simulations are able to attain, in order for locality and its consequences 
to become fully manifested. 
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APPENDIX A 



We give here more details and precise statements of our rigorous results. We first explain 
briefly the justification of the scaling relations in eq. Q that are employed in our argu- 
ment. A basic observation of experiments and numerical simulations^ ^^ * ^^ * ^ ^ on turbulent 
flow is the scaling of structure functions of velocity-increments Su{r; x) = u(x+r) — u(x) : 



(|5u(r)r> 



i/p 



Apr"" 



(A-1) 



for rjp <^ r L with rjp a short-distance viscous length-scale. As earlier, we shall 
interpret (■) as a space- average. We assume that the scaling range becomes longer as 
z/ ^ 0, so that rip 



0. A weaker form of "scaling" implied by (A-1) is then that 

log ||5u(r)|| 



lim inf 

r~>0 



logr 



(A-2) 



with 



I |p^i/p ^]2g Lp-norm. If we assume also that (|u|p) < oo, then relation 
(A-2) with < cTp < 1 implies that u is Besov regular with maximal Besov exponent cTp 

for all s < CTp but for no s > dp. E.g. geelSESESl fQj^ 



of order p. That is, u E B 
more background on Besov spaces and their relations to turbulent structure functions 



Assuming validity of (A-1) [or even (A-2)], a Paley-Littlewood criterion for Besov 
regularity implies a rigorous version of the second scaling relation in eq.(|4]): 



loe llul^'lll 

1- ■ r " \\p 

ap = lim mt 



(A-3) 



K/2<\k\<K 



\og{l/K) ' 

;*'*^'^u(k). See Triebel^^, Theorem 3.5.3(i). Here the norm on 



with u[^](x) = 

wavenumbers must be taken to be |k|oo = maxj/cj., ky, k^} or |k|i = + \ky\ + \kz\ and 
not the Euclidean norm |k|2 = (A;^ + ky + klY^"^. As we shall see below, other delicate 
aspects of Fourier analysis also rule out the use of the Euclidean norm for a rigorous 



proof of locality bounds. Note that (A-3) is actually more general than (A-2) and gives 



the maximal Besov exponent of u for any dp G M and p > 1. The first scaling relation in 



eq.(|4]) has a similar rigorous statement for Op > 0, 



. log||u>^||p 
(T„ = hmmf — — ; — , . 

P K^oc \og{l/K) ' 



(A-4) 



which follows from results on the strong summability of Fourier partial sums that char- 
acterize Besov spaces. See Triebel'^, section 3.73, Theorem l(i). Finally, u G i?p'°°(T^) 
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iff Vu G B'p-^'°°(T^) and then the result 



l = liminf '°fl^';;""' (A-5) 



for dp < 1 follows from (A-3). 

We now state and prove a rigorous analogue of inequality ([s]), implying UV locality 
of the spectral SGS flux: 

Proposition 1. If u has maximal Besov index a^p of order 3p forp > 1, then 

||n^(u,u,u[«i)llp<Cpiri"''Q~2'' 

when Q > AK, for any a < a^p. 

Proof of Proposition 1: The starting point of the argument is the identity 
derived in the text. Then by the Holder inequality 

We next exploit the Lp-boundedness of Fourier multipliers for polygonal partial summa- 
tion to obtain 

for some constant C > 0. See e.g., Krantzl^Sl^ Theorem 3.4.5. Note that such an estimate 
cannot be obtained with the Euclidean norm on wavenumbers, due to the failure of 
boundedness of Fourier multipliers for the ball. See Fefferman (1971)^3. However, we 
can approximate the sphere with any convex polygon with an arbitrary but fixed number 
of sides and our rigorous results will hold (see, for example, Krantzl^, and other references 
on the subjectSSESl). Using again the Holder inequality gives 

||n;,(u,U,u[Q])||, < C||Vu<^||3,||u[?-^-«+^]||3,||u[«l||3,. 



Writing u^z -^'Q+^l = u^? ^ — w>^^^ ^ we can use (A-4) to obtain the bound 
||u[?-^'«+^]||3p < ||u>^-^||3p + ||u>«+^^l3p < C'p [{Q - 2Kr^ + {Q + K)-^] < C^Q-^ 
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for any a < a-^p and Q > AK. Using the similar estimates for a < a^p 



from (A-5) and (A-3) yields the final inequality. 



□ 



Finally, we state and prove a rigorous analogue of inequality (13), implying IR locality 
of the spectral SGS flux: 

Proposition 2. If u has maximal Besov indices a2 and of order 2 and 4, resp., and 



if the strong scaling relation (12) holds for the energy spectrum, then 



||nK(u,u[^U)l|i < ciri/2-'-pi/2-' 

when P < K/2, for any a < a2 and any a' < (T4. 
Proof of Proposition 2: We begin with the identity 



<K 



(-1 



[P] [K-P,K+P]\<K [P] / [K-P,K+P\\<K 



U 



equivalent to ( 11 ) in the text. By applying the Holder inequality and the Lp-boundedness 



of Fourier multipliers for polygonal partial summation, then, similarly as above, we obtain 

||n^'(u,u[^l,u)||l < C||Vu<^||4||u[^]||4||u[^-^'^'+^]||2. 

Since Hu'-'^"^''^"'"^] ^ = ||u^^^~'^||2 — ||u^-^^+^||2, we obtain from the first relation in Q 

||^[A'-p,A-+p]||2 = c[{K~ P)-^''^ -{K + P)-2-^] < C'PK-^'^"'^ 
for any P < K/2 and any a < (T2. Finally, employing the estimates for a' < (T4 

||Vu<^||4 < C'K'-'^', ||u[«l||4 < CQ-""' 



from (A-5) and (A-3) with p = A yields the final inequality. 



□ 



APPENDIX B 



We provide here some details of the derivation of equations (23) and (24) in the text. 



Thus, we take k = K — qv and p = K + qu and evaluate all expressions to leading order 
in the small parameter 5 = q/ K ^ 1. 
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We begin by evaluating the interaction coefficient Bkqp in (21). Using the law of 
cosines, k"^ = + — 2pqcosa, one obtains 



X = cos a = 



+ — k"^ 
2pq 



zpq 

K\l + 6u)'' + K^^ - K\l - 6v)'' 

2- K{l + 6u)- K5 
= u + v + 0{5) 

In the same manner, 

y = cos P = —{u + v) + 0(6) and z = cos'j = 1 + 0(5'^). 
Thus, up to terms of relative order 0{6), 

Bkqp = ATC^kYqixy + z'') ~ An^K^q{l - w^) 

with w = u + V. 

We next substitute 



U{p; t\s) = U{K- t\s) + qu^U{p; t\s) \^^^ + 0{q^ 



G{p- t\s) = G{K- t\s) + qu^G{p- t\s\^^ + 0{q^) 

and the similar expansions for U {k; t\s), G{k; t\s) into the time-integral in (21 ). Employ- 
ing also the inert ial-range scaling relations (22) and the anti-symmetry of the integrand 
in k and p, one obtains up to terms of relative ^---i"- /^f^^ ■ 



tive order 0{6) : 

[ ds[G{k;t\s)U{p;t\s) - G{p;t\s)U{k;t\s)]U{q;t\s) = q{u + v) 
ho 

X f ds{gie"'K^'\t - s))^\—e^"p-'"hie"^r?l^ 
Jto I ^dp^27t 

G 



. G^ „ 



To obtain the last line the integration variable was changed to r = e^^^k'^^^{t — s) and 
the resulting integral evaluated to leading order to be / = y/i + |/2, with 

/•OO /"OO 

^(r)r(r) dr = 0.76, h = [^'(r)r(r) - ^(r)r'(r)] r dr = -0.19. 
Jo Jo 
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These latter two integrals were numerically evaluated by Kraichnan'^, eq.(6.2). 
Finally, one obtains from (21): 

T{k,q,p) ^ -ew{l -w^)C^IK~^/'^q~^/\ 

with I = 2.66 and C = 1.77 the Kolmogorov constant from the ALHDIA closure (see 
Kraichnanfl^, eq.(6.1)). This is exactly equation (23) in the text. 

We now substitute this result into (20) to obtain (24). The integral over u and v is 
easily evaluated by the change of variables w = u + v, t = {u — v)/2, giving 

1 pl — v pi nw/2 pi 2 

dv / duw{l — w^)= / dw / dtw{l — w'^)= / dw w'^ {1 ~ w"^) = — . 
Jo Jo J-w/2 Jo 15 

It may be useful to discuss briefly the relationship of our heuristic derivation of the 
4/3-type scaling of Uk based on decorrelation effects in Section II.C.l and the 

derivation using Kraichnan's ALHDIA closure presented in this appendix. The separate 
input and output terms in the ALHDIA expression for T{k,q,p) each yield, on order of 
magnitude, 

n^/°"*(u,u[QU) ^ k'q't{k)e{k)e{q) 

where t(K) ~ e^^^^K^"^^^ is the eddy-turnover-time at wavenumber K. This factor arises 
from the time-integration in eq.(21) and represents the effect of decorrelation in Kraich- 
nan's closure. Using the relations that KE{K) ^ [^^(fx)]^ and Kt{K) ^ l/5u{^K)■, this 
result becomes 

However, because of the near cancellation between input and output terms for Q <^ K, 
the final ALHDIA result is smaller by a factor of Q/ K = ix/iQ, giving 



H™^(u,u[«U) ^ KQ'r{K)E{K)E{Q) ^ 

This corresponds very closely to our result and derivation in Section II.C.L The scaling 
of the partial flux H7^(u, u^^\ u) in our eq.(14) there is due to the difference, 

where the two terms scale similarly to each other and nearly cancel. This cancellation 
yields an additional factor of ix/^Q, so that our final result scales as 
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